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Section  1 
INTRODUCTION 


1. 1 DESCRIPTION  IF  GASFORMING  PYROTECHNIC  REACTIONS 


During  the  burn  of  a pyrotechnic  composition,  the  solid  is  transformed  into  a gas  by 
pyrolysis,  chemical  reaction  and  evaporation.  Even  if  the  transformation  occurs 
through  the  formation  of  an  intermediate  liquid  phase,  there  exists  a certain  position, 
called  the  burning  surface,  where  the  gas  is  being  formed.  The  mechanism  of  the 
solid-phase  decomposition  at  the  burning  surface  is  considered  to  be  the  rate- 
controlling step,  the  gas  being  formed  at  the  rate 


pr  = pB  exp  ( -E/RT  ) 

s 


(1) 


where 

r = burning  rate  (cm/sec) 

B = the  pre-exponential  frequency  factor 
E = the  activation  energy 
R - the  gas  constant 

T = the  temperature  at  the  burning  surface 

s 3 

p = the  density  (gm/cm  ) 

Reactions  occurring  in  the  solid  phase  would  have  changed  the  temperature  at  the 
burning  surface  to  Tf , but  actually  a temperature  T greater  than  T'  is  obtained 

S 8 3 

there  because  heat  is  conducted  back  to  the  solid  from  the  gas.  The  essential  theore- 
tical problem  is  to  discover  how  T depends  on  the  gas-phase  process. 

o 
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Tracer  mixtures  constitute  not  homogeneous  materials  but  rather  composites  that 
contain  both  discrete  fuel  and  oxidizer  particles,  the  decomposition  of  which  is 
governed  by  Eq.  (1).  For  the  steady-state  burning  of  a composite,  it  is  reasonable 
to  assume  that  the  linear  rates  of  regression  of  the  oxidizer  and  the  fuel  surfaces  are 
about  equal: 


"l  rt  - "o  r„ 

"t  Bf  <-VRT„f>  ” "f  B0  e*P  <-VRTso> 


(2) 


where  the  subscripts  f and  o refer  to  the  fuel  and  oxidizer,  respectively. 


In  general,  the  preexponential  factor  and  the  activation  energy  of  the  fuel  and  oxidizer 
are  not  equal.  It  follows  that  Tgf  and  T must  also  differ.  Because  of  the  difference 
in  pyrolysis  rates,  one  constituent  gasifies  relatively  faster,  and  the  slower  burning 
component  is  left  protruding  from  the  surface.  The  faster  pyrolysing  constituent  is 
therefore  in  contact  with  a cooler  region  of  the  flame  than  is  the  slower  burning  con- 
stituent. Eventually,  particles  of  the  slower  burning  constituent  may  be  ejected 
entirely  from  the  surface  and  be  entrained  in  the  hot  gas  stream  where  they  will  under- 
go further  reaction. 


1.  2 DESCRIPTION  OF  TRACER  REACTION 


Tracer  reactions  are  well  described  by  the  above  model.  A recent  study  (Ref.  1)  of  a 
motion-picture  record  of  a burning  tracer  composition  has  shown  that  most  fuel  parti- 
cles (atomized  magnesium)  are  ejected  as  molten  globules  from  the  surface.  Because 
of  the  high  gas  temperature,  the  fuel  particles  vaporize  so  as  to  give  rise  by  subsequent 


oxidation  to  a hot  reaction  zone.  This  flame  front  is  narrow  because  of  the  high  rate 
of  the  gas-phase  reaction  as  compared  with  the  preceding  solid-phase  decomposition 
and  evaporation  steps.  The  temperature  and  the  proximity  of  the  flame  to  the  surface 
depends  (at  constant  ambient  pressure)  on  the  particle  size  and  the  stoichiometry. 

Tracer  compositions  are  usually  fuel-rich,  resulting  (as  will  be  shown  in  a later  sec- 
tion) in  an  equilibrium  gas  composition  containing  large  quantities  of  magnesium  vapor. 
This  combustible  gas  is  largely  responsible  for  the  luminous  plume  of  the  tracer 
bullet. 

In  describing  the  burning  of  magnesium  particles,  two  reaction  models  have  been  con- 
sidered (Ref.  2),  which  differ  mainly  in  the  treatment  of  the  condensed  oxide  formed.  Ex- 
perimental evidence  shows  ihat  some  magnesium  particles  will  react  with  an  oxidizer 
that  surrounds  them  with  a crust  of  oxide.  The  reaction  rate  is  determined  by  the  dif- 
fusion through  this  crust  which,  owing  to  the  higher  density  of  the  oxide  relative  to  the 
parent  metal,  is  porous  and  readily  permits  the  counterdiffusion  of  gaseous  oxidizer 
and  fuel.  The  end  product  is  an  empty  bubble  that  is  usually  a part  of  a cake  which 
which  was  formed  in  situ  by  the  interaction  with  the  oxide  growth  from  the  neighboring 
particles.  These  oxide  bubbles  are  shown  in  Fig.  la. 

Most  of  the  magnesium  becomes  entrained  in  the  stream  of  gaseous  oxidizer.  While 
the  rate  of  heat  input  in  the  flame  zone  suffices  to  cause  the  magnesium  to  evaporate, 
the  resulting  vapor  oxidizes  only  partially,  because  of  the  fuel  rich  nature  of  the  mix- 
ture, so  that  the  oxide  particles  remain  detached  from  the  metal  surface.  These  oxide 
particles  have  colloidal  dimensions  and  are  entrained  in  the  gas  stream  with  the  other 
products  of  the  reaction.  The  reaction  rate  is  controlled  by  the  vapor  phase  diffusion 
of  magnesium  to  the  reaction  site.  This  reaction  model  best  describes  the  process 
which  we  have  observed  and  is  illustrated  in  Fig.  lb. 

The  the°ret'cal  problem  in  describing  tracer  reactions  is  to  learn  the  effects  cf  * ue 
stoichiometry  and  particle  size  on  the  ratu  of  energy  release  in  the  fU;me  and  hence 
the  energy  partition  between  the  heat  deposited  into  the  solid  and  the  heat  entrained 
in  the  exhaust  gas.  The  effect  of  reducing  the  fuel  particle  size  (at  constant  composition) 
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120  fi  I 1 85x 

Fig.  la  Residual  Magnesium  Oxide  Bubbles.  The  product  was  obtained  by 

burning  a tracer  bullet  in  an  argon -filled  chamber  and  by  separating 
the  granular  solids  from  the  colloidal  dust  which  adhered  to  the 
reactor  walls.  These  hollow  bubbles  remain  after  the  magnesium 
particles  are  coated  with  an  oxide  crust  followed  Dy  flash  evapora- 
tion of  the  metal  through  a vent  hole.  The  rate  of  this  process  is 
too  slow  to  account  for  the  tracer  reaction. 

HIGH 


Fig.  lb  Predominant  Burning  Mechanism  of  Magnesium  Farticles  (Ref,  2) 
Fig.  1 Models  for  Burning  of  Magnesium  Particle 
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will  he  to  narrow  the  flame  zone  and  to  increase  the  energy  partition  and  hence  the 
regression  rate  of  the  surface.  The  result  will  be  a brighter  plume  of  shorter  dura- 
tion. Reduction  of  the  fuel  (at  constant  particle  size)  may  reduce  thermal  conductance 
in  the  solid,  reduce  the  regression  rate,  and  prolong  but  dim  the  trace.  At  the 
extremes  in  composition,  the  trace  quality  will  become  insensitive  to  fuel 
characteristics. 


1.3  OBJECTIVE 


The  purpose  of  this  study  is  to  construct  an  analytical  model  that  will  put  the  preceding 
quantitative  description  on  a quantitative  basis.  An  understanding  of  all  the  factors 
that  determine  the  energy  partition  is  beyond  the  state  of  the  art,  and  for  this  reason 
the  baseline  conditions  were  to  be  determined  by  experimentally  determined  regres- 
sion rates.  Inasmuch  as  this  study  is  of  preliminary  but  otherwise  perfectly  general, 
nature,  the  effect  cf  binders  and  of  color  intensifiers  on  tracer  performance  has  been 
omitted  from  the  theoretical  analysis.  However,  methods  by  which  such  variables 
can  be  accounted  for  and  the  effect  of  such  additives  on  the  equilibrium  composition 
of  the  product  grades  will  be  described  in  this  report. 


1.4  TECHNICAL  APPROACH 


The  rate  constants  of  the  decomposition  of  the  components  were  determined  experi- 
mentally. The  thermochemical  properties  were  collected  from  published  sources. 
The  energy  partitions  were  calculated  from  several  compositions  of  a single  particle 
size  of  the  magnesium  using  experimentally  determined  burning  rates  of  spinning 
7. 62-mm  tracer  bullets.  Thermal  conductivities  were  measure  i using  the  thermal 
diffusion  method  (Ref.  3). 
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The  analytical  problem  concerns  the  heat  transfer  into  the  endothermally  decomposing 
solid  with  moving  boundaries.  The  gas -phase  process  was  treated  by  examining  the 
experimental  evidence  for  the  case  of  fuel  vapor  diffusing  through  an  oxide  shell  and 
for  the  case  of  a detached  reaction  zone  surrounding  the  droplets.  It  was  found  that 
only  the  detached  reaction  zone  model  (Fig.  lb)  could  account  for  the  experimentally 
observed  reaction  rates . 
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Section  2 

EXPERIMENTAL  STUDY 

2.1  DE  COM  POSITION  KINETICS  OF  STRONTIUM  NITRATE 

2. 1. 1 Literature  Survey 

The  kinetics  of  the  decomposition  of  inorganic  nitrates  has  never  been  studied  exhaus- 
tively, particularly  with  respect  to  the  alkaline  earths,  and  recent  studies  are  almost 
completely  lacking.  This  situation  is  surprising  in  view  of  the  wide  range  of  applica- 
tions of  inorganic  nitrates  in  the  pyrotechnics  and  fertilizer  industry. 

The  question  that  had  to  be  answered  in  this  study  was  whether  the  nitrate  or  the 
nitrite  was  the  more  stable,  i.e. , whether  the  decomposition  of  the  nitrate  could  be 
treated  as  a first-order  reaction  with  respect  to  its  decomposition  or  whether  some 
other  decomposition  step  was  rate  determining.  Although  no  rate  data  have  ever  been 
published,  the  equilibria  have  been  studied  by  Centnerzwer  (Ref.  4),  who  found  that 
the  nitrite  decomposes  slowly  below  its  melting  point  of  264°C: 

380°C 

Sr(NO„)  — SrO  + N„0  (3) 

1 atm 

The  nitrous  anhydride  disproportions. tes  in  a fast  step: 

N203=N02  + NO  (4) 


so  that  the  resulting  dioxide  can  react  with  more  nitrite  in  another  fast  reaction: 

Sr(N02)2  + 2 N02=SSr(N03)2  + 2 NO  (5) 
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Because  the  rate  of  oxidation  of  the  nitrite  is  greater  than  the  rate  of  its  dissociation, 
the  nitrate  will  be  the  more  stable  specie.  This  was  shown  (Ref.  5)  by  the  stability  of 
the  nitrate  up  to  580°C.  This  stability  of  the  nitrate  was  confirmed  by  differential 
thermal  analysis  (6).  Findings  that  the  nitrate  decomposes  in  several  steps  starting 
at  280°C  (7)  must  have  been  due  to  test  samples  that  were  contaminated  by  the  nitrite. 


Alkali  nitrates,  by  contrast,  form  an  intermediate  nitrite  on  decomposition  (Ref.  G)  1 

I 

which  was  explained  by  the  fact  that  alkali  nitrates  dissociate  above  the  melting  point  I 

and  hence  have  two  degrees  of  freedom,  whereas  the  alkaline  earth  nitrates  diasociale  j 

below  the  melting  point  (SrfNO^)^,  618°C,  Ref.  6)  and  therefore  have  only  one  degree  j 

of  freedom  (Ref.  7).  < 

We  conclude  from  this  study  that  strontium  nitrate  is  stable  with  respect  to  its  decom-  | 

position  into  the  nitrite  and  that  the  decomposition  is  most  likely  a first  order  reaction  \ 

with  respect  to  the  nitrate:  j 

* 

2 Sr(NOa)2 - 2 SrO  + 4 NO  + 3 02  (6)  J 

i 

* 

2.1.2  Experimental  Study  1 


J 

Strontium  nitrate  was  decomposed  isothermally  in  a thermal  gravimetric  analyzer  \ 

(TGA).  It  was  found  that  unless  air  (containing  C02)  was  excluded,  the  oxide  product 

was  converted  into  the  carbonate.  Experimental  data  became  reproducible  when  the 

decomposition  was  performed  in  helium.  The  resulting  weight  losses  equivalent  to 

a s 

Sr(NO„)„  were  used  to  compute  the  ratio where  a was  the  original  weight  of  the  ! 

w fc  3 " X | 

nitrate  at  time  zero,  and  x was  the  weight  loss.  Figure  2 shows  these  normalized 
weight  loss  data  plotted  as  function  of  time  for  five  temperatures,  the  slope  of  the  line 
being  the  first-order  rate  constant  k^,.  The  data  were  trected  by  a Tymshare  least 
squares  routine.  The  results  are  shown  in  Table  1. 
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Fig.  2 Decomposition  Rate  Data  for  Strontium  Nitrate 
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Table  1 


i 


I 

s 

t 

5 


RESULTS  OF  LEAST  SQUARE  ANALYSE  OF  STRONTIUM  NITRATE 
DECOMPOSITION  RATE  DATA 


Temperature 

(°K) 

Computed 

i n ® . - it  t 

(min  ) 

Std.  Dev. 

■ a 

ln  a - x “t* 

in  in  ■■ 
a - x 

917 

0.00470  + 0.01316  t 

0.0132 

0.0151 

929 

- 0.01885  + 0.02476  t 

0.0248 

0.0276 

943 

0.00582  + 0.03613  t 

0. 0361 

0.0221 

960 

0.05458  + 0.05363  t 

0. 0536 

0.0454 

973 

- 0.02179  + 0.08559  t 

0.0856 

0.0216 

The  ra<e  constants  k^.  were  correlated  with  the  reciprocal  temperature  as  shown  in 
Fig.  o . The  slope  of  the  line  is  proportional  to  the  energy  of  activation,  E . We  con- 
clude that  the  linear  dependence  of  the  decomposition  data  on  time  justifies  the  assump- 
tion of  a first-order  rate  law  which  is  expressed  as  follows: 

C = Cf'*  (7) 

and 

Ry  = Ze  RT  (8) 


From  Fig.  2 


E = 55,400  ± 6000  cal/mole  and  Z = 1 x io11  min  * 

2.2  PARTICLE  SIZE  DISTRIBUTION  OF  MAGNESIUM 

Figure  4 shows  u scanning  electron  micrograph  of  the  Frankford  Arsenal  magnesium. 
The  size  distribution  was  determined  by  particle  count,  and  the  characteristic  size  was 
found  to  be  120  p. 
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2. 3 THERMAL  DEGRADATION  OF  POLYVINYL  CHLORIDE 


Unlike  the  dissociation  of  strontium  nitrate,  the  thermal  response  of  PVC  is  complex. 
Figure  5 shows  the  typical  appearance  of  low -temperature  isothermal  TGA  data.  The 
sigmoid  shape  of  the  curve  is  typical  of  chain  reactions,  which  for  PVC  may  be  vis- 
ualized as  follows.  The  carbon -chlorine  bond  ie  the  weakest  bond  in  the  polymer  so 
that  it  is  first  to  rupture.  This  is  followed  immediately  by  an  abstraction  of  a hydro- 
gen atom  from  an  adjacent  carbon  and  the  formation  of  a double  bond  in  the  chain. 

HHHHH  HHHHH 
I l I l l I i t I i 

-C-C-C-C-C *~C  - C = C - C - C ~ 4 HC1  (9) 

I I i I i I it 

Cl  H Cl  H Cl  H Cl  H 

The  chlorine  atoms  in  the  0 position  to  the  double  bond  are  now  more  likely  to  break 
away  from  the  carbon  atoms  to  which  they  are  attached  than  the  other  chlorine  atoms. 
Hence,  once  one  chlorine  is  removed  (requiring  a relatively  high  activation  energy, 
and  consequently,  a lengthy  induction  period),  the  others  come  off  easily,  leaving 
behind  a polyunsaturated  structure  which  is  prone  to  further  reaction.  It  may  cross- 
link, rupture,  or  form  cyclic  compounds.  This  explains  why  the  analysis  of  the  vol- 
atile decomposition  products  show  in  addition  to  HCI  also  ethylene,  propylene  ar.d 
benzene  (Ref.  8).  The  remainder  is  a polyacetylene  char  which  is  subject  to  further 
attack  by  oxidizers. 

The  above  steps  (except  for  the  initial  activation)  are  strongly  exothermic,  and  it  has 
been  found  that  when  PVC  is  heated  in  a confined  volume,  it  will  decompose  explo- 
sively (Ref.  10).  Using  differential  scanning  calorimetry,  a value  of  -72  -tlO  cal/g 
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VO 


Fraction  of  HC1  Remaining  in  Polyvinyl  Chloride  at  252  *C  as  a Function  of  Time 


was  found  for  the  heat  of  decomposition  of  PVC,  This  value  is  small  compared  with 
otter  thermal  quantities  that  control  the  reaction.  It  has  therefore  been  ignored  in 
this  study. 

Dehydrochlorination  of  PVC  appears  to  follow  a 3/2  order  reaction  rate  law  (Ref.  9). 
The  equation  for  such  a reaction  is 


= k C^2 
dt 

which,  upon  integration,  yields 

1 1 


(10) 


where  C is  the  fraction  of  HC*  remaining  in  the  sample  at  time  t,  and  k^,  is  the  specific 

reaction  rate  constant.  The  fit  of  the  experimental  data  as  obtained  from  isothermal 

TGA  runs  is  shown  in  Fig.  6.  An  Arrhenius  plot  based  on  the  resulting  rate  constants 

(see  Fig.  7)  yielded  an  activation  energy  of  10,750  cal/ mole  and  a pre -exponential 
8 -1 

factor  of  10  min  , using  Eq.  (8). 


Table  2 

HATE  CONSTANTS  FOR  THE  THERMAL  DEGRADATION 
OF  POLYVINYL  CHLORIDE 


kT 

Temperature 

C‘l/2/min 

fC> 

252 

0.11 

272 

0.14 

294 

0.21 

We  conclude  that  the  presence  erf  PVC  in  the  tracer  mixture  will  complicate  the  con- 
struction of  a tracer  reaction  model  in  that  the  presence  of  two  decomposing  solids 
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Thermal  Degradation  of  Polyvinyl  Chloride  Test  of  3/2  Order  Reaction  Rate  Law 


Fig.  7 Thermal  Decomposition  of  Polyvinyl  Chloride  — Plot  of  Reaction 
Rate  Constants 


a 


will  require  a two-zone  decomposition  model;  moreover,  this  model  is  characterized 
by  an  exothermic  low-temperature  reaction  followed  by  a high -temperature  endother- 
mic reaction.  Such  models,  when  programmed  for  computer  analysis,  tend  to  be 
costly  in  terms  of  machine  time  because  of  multiple  zones.  To  simplify  this  analysis, 
tills  study  will  be  confined  to  a two-con^  jnent  system  - magnesium  in  strontium  ni- 
trate with  weight  fraction  and  particle  size  of  the  magnesium  as  variables. 

2.4  EQUILIBRIUM  CALCULATIONS  FOR  VARIOUS  TRACER  COMPOSITIONS 

Several  tracer  compositions,  both  fuel-rich  and  stoichiometric  (Table  3). were  subjected 
tc  a computer  study  of  the  reaction  products.  The  results  are  reported  in  terms  of 
flame  temperature,  mole  fraction  In  the  gas  phase,  mol  fraction  in  the  solid  phase , 
mole  fraction  in  the  solid  phase  based  on  the  total  composition,  and  the  percentage 
of  product  moles  to  total  moles  of  elements. 

It  will  be  noted  from  the  results  shown  in  Tables  4 through  9 that  the  product  is 
chiefly  composed  of  CO,  H2>  gaseous  Mg,  Ng,  and  SrCl2  and  a solid  phase  which 
consists  chiefly  of  MgO.  The  addition  of  PVC  tends  to  shift  the  equilibrium  from  SrO 
to  SrClgl  in  contrast,  increasing  tb  j amount  of  oxidizer  increases  the  water  vapor  at 
the  expense  of  the  hydrogen.  In  fuel-rich  systems,  unreacted  magnesium  and  hydro- 
gen are  the  chief  reactive  components  if  it  can  be  assumed  that  the  available  o:Jdizer 
is  effectively  depleted  within  the  bullet  cavity.  This  finding  would  have  the  important 
simplifying  implication  (for  fuel  rich  systems)  that  the  critical  parameters  affecting 
the  performance  of  the  tracer  munition  are  the  kinetic  and  thermal  transport  processes 
occurring  in  the  solid  phase. 
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Table  3 

PRODUCT  DESCRIPTION  FOR  VARIOUS  TRACER  FIXTURES 


Tracer  Composition 
PVC/Mg/Sr<NO,)_ 
(Wt  %)  32 

Mixture 

Description 

Flame 

Temperature 

no 

Moles  Comb.  Prod  _ 
Moles  of  Elements  x 

m 

2000 

28.0  - 9.1  - 63.0 

Stoichiometric 

2980 

43.66 

42.01 

41.92 

19.8  - 12.7  - 67.5 

Stoichiometric 

3060 

44.32 

41.45 

41.12 

30.0  - 19.4  - 50.0 

Fuel  rich 

2840 

49.44 

49.18 

43.99 

22.1  - 28.5  - 49.6 

Fuel  rich 

2980 

51.84 

51.92 

51.95 

17.0  - 28.0  - 55.0 

R-284 

3070 

50.07 

50.02 

17.0  - 33.0  - 50.0 

Fuel  rich 

3070 

53.24 

53.37 

. — - - 

1 
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Table  4 

EQUILIBRIUM  COMPOSITION  OF  REACTION  PRODUCTS 
TRACER  COMPOSITION  (PVC  - Mg  - Sr  (NO^  ) 28.0  - 9.1  - 63.0wt% 


Table  5 


EQUILIBRIUM  COMPOSITION  OF  REACTION  PRODUCTS 
TRACER  COMPOSITION  (pVC  - Mg  - Sr  (N03)2  ) 19.8  - 12.7  - 67.5  wt% 


Mole  Fraction 

Product 

3000"  K 

2500°  K 

2000°  K 

Gas 

Solid 

Gas 

Solid 

Gas 

Solid 

CO 

0. 1953 

0. 1263 

0. 1088 

C02 

0. 1330 

0. 2540 

0, 2763 

Cl 

0. 0015 

0. 0003 

— 

HCi 

0. 0025 

0. 0023 

0.  0015 

Mg  Cl 

0.  0003 

— 

— 

H Mg  O 

0. 0148 

0. 0007 

— 

HO 

0. 0461 

0. 0087 

0.  003 

«2 

0. 0317 

0.0187 

0.  0205 

h2° 

0. 1555 

0. 2290 

0.  2355 

Mg 

0.  0357 

0. 0003 

— 

NO 

0.0119 

0. 0021 

— 

N2 

0. 1896 

0.  2256 

0.  2294 

O 

0. 0251 

0. 0011 

— 

°2 

0. 0497 

0. 0053 

— 

SrC^ 

0. 1073 

0. 1256 

0. 1277 

MgO 

0. 7542 

0.7862 

0.  7876 

(Fraction  of  Total 
Mole  Product) 

(0. 1995) 

(0.2526) 

(0. 2554) 

jSrO 

0.  2458 

0.2138 

0.  2124 

(Fraction  of  Total 
Mole  Product) 

(0.065) 



(0.06869) 

(0.  0685) 
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Table  7 


EQUILIBRIUM  COMPOSITION  OF  REACTION  PRODUCTS 
TRACER  COMPOSITION  (PVC  - Mg  - Sr  (NOgIg  j 22.1  - 23.5 


- 49.6  wt  % 


3000* K 


Mole  Fraction 
2500“  K 


2000*  K 


0. 2638 
0.0219 
0.  0005 
0.  0015 
0.  0009 
0.  0272 
0.0103 
0. 1070 
0.0640 
0. 2928 
0. 0011 
0.  1116 
0. 0031 
0. 0007 
0.0938 


MgO 

(Fraction  of  Total 
Mole  Product) 


Solid 


Gas 

C-Ifj  Ga3 

0. 3237 

0. 3269 

0. 0010 

— 

0. 0003 

— 

0. 0002 

— 

0. 0023 

— 

0.  2110 

0.2178 

0.  0041 

— 

0.  2217 

0.2179 

0. 1275 

0.1283 

Sn; 

! (Fraction  of  Total 
I Mole  Product) 


(0. 01438) 


(0.01337) 


Tabic  8 


EQUILIBRIUM  COMPOSITION  OF  REACTION  PRODUCTS 
TRACER  COMPOSITION  (p VC  - Mg  - Sr  (NO^  ) 1?  - 28  - 55wt% 


Mole  Fraction 

Product 

3000*  K 

2500* K 

2000aK 

Gas 

Solid 

Gas 

Solid 

Gas 

Solid 

CO 

0. 2631 

0. 3379 

0. 3434 

O 

o 

0. 0282 

0.  0019 

— 

a 

0.0006 



h a 

0.  0016 

0.  0004 

— 

Mg  a 

0.0008 

0.0002 

— 

H Mg  o 

0. 0260 

0. 9024 

— - 

HO 

0.  0128 

— 

— 

H2 

0. 0980 

0. 2175 

0. 2287 

v 

0.  0758 

0.0073 

— - 

Mg 

0.  2265 

0. 1271 

0. 1190 

NO 

0. 0017 

— 

— 

N2 

0. 1640 

0. 1922 

0, 1943 

o 

0.  0039 

— 

— 

°2 

0.  0012 

— 

— 

SrCl2 

0.0956 

0. 1130 

0. 1145 

MgO 

0. 8732 

0.9011 

0. 9028 

(Fraction  of  Total 
Mole  Product) 

(0. 3088) 

(0.4011) 

(0.4075) 

SrO 

0. 1268 

0.0989 

0. 0972 

(Fraction  of  Total 
Mole  Product) 

(0. 0448) 

(0. 04405) 

(0. 0439) 

Table  9 


EQUILIBRIUM  COMPOSITION  OF  REACTION  PRODUCTS 
TRACER  COMPOSITION  (pVC  - Mg  - Sr  (NOg)2  ) 17  - 33  - 50wt% 


Mole  Fraction 

Product 

3000* K 

2500“  K 

2000“  K 

Gas 

Solid 

Gas 

Solid 

Gas 

Solid 

CO 

0. 2352 

0. 2785 

0. 2804 

C°2 

0. 0165 

0. 0007 

— 

Cl 

0.0004 

— 

— 

H Cl 

0. 0013 

0.  0002 

— 

Mg  Cl 

0.0009 

0.0002 

— 

H Mg  0 

0.  0262 

0.  0022 

— 

HO 

0.0084 

— 

— 

H2 

0.  0995 

0. 1823 

0. 1869 

H2° 

0.  0505 

0.0026 

— 

Mg 

0.  3453 

0. 2969 

0. 2949 

NO 

0.0010 

— 

--- 

N2 

0. 1290 

0. 1436 

0. 1442 

O 

0.0026 

— 

~ 

^2 

0.  0005 

— 

^ _ 

Sr  Clg 

0.  0826 

0.0928 

0.0935 

MgO 

0.8880 

0.9120 

■■ 

0.9  IS? 

(Fraction  of  Total 
Mole  Product) 

(0. 2620) 

(0.  3335) 

(0. 3360) 

SrO 

0. 1120 

0. 0880 

0. 08684 

(Fraction  of  Total 
Mole  Product) 

(0. 0331) 
— ■ 

(0. 03219) 

. _ . 

(0.  0320) 
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2.5  THERMAL  TRANSPORT  PROPERTIES  OF  TRACER  MIXTURES 


The  experimentally  determined  thermal  transport  properties  are  listed  in  Table  10 
for  mixtures  of  magnesium  with  strontium  nitrate.  These  data  are  plotted  in 
Figure  8.  It  will  be  observed  that  the  conductivities  of  magnesium -rich  mixtures 
are  significantly  below  those  of  bull;  magnesium  (k  * 0. 38  cal/cm  sec  ®K). 


Table  10 

THERMAL  TRANSPORT  PROPERTIES  OF  TRACER  MIXTURES 
Determined  By  Thermal  Diffusivity  Method  (Ref.  3) 


Heat  Capacity 
(cal/’K  g) 


Thermal  Conduc- 
tivity 

{ca2/cm~sec)®K 


1.57  x 10~3 
1.83  x 10~3 


THERMAL  CONDUCTIVITY  (CAl/CM-8IC*K) 


Section  3 

ANALYTICAL  STUDY 


3. 1 SOLID-PHASE  PROCESS 

The  mathematical  representation  of  the  thermal  transport  into  the  solid  is  pictured 
below. 


The  location  of  the  surface  is  denoted  by  s(t).  It  regresses  at  a constant  speed  , 

steady-state  conditions  being  assumed. 

The  thermal  transport  process  consists  of  the  heat  fine  (H)  into  the  surface  and 
heat  absorption  by  the  thermal  decomposition  of  Sr(NOg)2.  The  fusion  of  mag- 
nesium and  the  presence  of  the  polyvinyl  chloride  are  ignored.  To  aid  in  the  visualiza- 
tion of  the  thermal  processes,  the  individual  steps  are  shown  separately  in  Table  11. 
The  thermal  transport  in  the  solid  is  described  by 

CP  * If  * * ~T+  pQ  t>o  on 


where 


= FZe 


E 

RT 


DF 

at 
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(12) 


and 


H + it 


dT 

dx 


0 at  x = s (t) 


03) 


Hie  boundary  conditions  are 


T <s  ft) ) = Ts 
F H) ) s o 


S (t) 


(14) 


Q is  the  heat  of  dissociation  of  strontium  nitrate.  The  kinetic  parameters  Z and  £ are 
determined  experimentally.  The  fraction  F of  undecomposed  material  has  a fixed  value 
at  a given  distance  from  the  surface.  The  temperature  gradient  at  the  surface  is  deter- 
mined by  the  heat  flux  H and  by  the  thermal  conductivity  of  the  material  |Eq.  (13)J. 
These  equations  are  solved  numerically  by  a computer  program  that  gives  the  regres- 
sion rate  when  the  heat  flux  is  specified.  This  computer  program  is  listed  in  the  Appen- 
dix . A temperature  profile  is  associated  with  each  heat  flux. 


For  a given  heat  flux,  fuel-rich  mixtures  are  expected  to  yield  higher  regression  rates 
for  two  reasons:  (1)  the  endothermic  heat  of  decomposition  varies  linearly  with  oxidizer 
content;  (2)  a higher  metal  fraction  brings  about  higher  thermal  conductivity.  Experi- 
mental data  on  regression  rates  as  reported  in  Table  12  are  in  agreement  with  this 
trend.  The  retarding  effect  of  adding  a binder  (PVC)  to  the  fuel-rich  mixture  is  evident 
from  the  data  on  the  R-284  mixture. 
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Table  12 


EXPERIMENTALLY  DETERMINED  REGRESSION  RATES 
(furnished  by  Frankford  Arsenal) 


Burning  Rate 

Description  of  Composition  Re^e'8‘“"  «atf  (om/sec) 

n = Spin  Rate  (rpm)  at  O = 40, 000  rpm 


Oxidizer  Rich  Composition 
28 .8%  Mg  71.2%Sr(NOs)2 

Stoichiometric  Composition 
36. 3%  Mg  63.  7%  Sr(NO,)0 

U Ci 

Fuel-Rich  Composition 
42.  8%  Mg  57. 2%  Sr(NOg)2 


R - 284 

55%  Sr(NO  ) 28%  Mg  17%  PVC 

3 * 


0.15  + 4.66  x 10  0 


0.31  + 3.91  x 10  fi 


0.41  + 3.71  x 10  D 


3.  2 GAS-PHASE  REACTION 


3. 2. 1 Introduction 


The  rate  of  energy  release  (heat  flux)  depends  on  the  gas-phase  reaction  which,  in  turn, 
is  governed  by  the  particle  size  of  the  fuel  and  the  gasification  mechanism.  The  overall 
tracer  burning  rate  depends  on  the  energy  partitioning  between  the  energy  deposited  into 
the  surface  and  the  total  energy  output  of  the  reaction.  Some  of  the  fuel-rich  mixture 
will  burn  outside  the  bullet,  but  this  is  of  no  interest  here.  The  enerpv  partitions 
were  not  determined  a priori,  but  were  obtained  from  the  experimentally  determined 
regression  rate. 


Inspection  of  a motion  picture  (Ref.  1)  of  burning  R-284  tr  acer  mixture  showed  that 
magnesium  particles  are  ejected  from  the  pyrolyzing  surface,  and  are  then  entrained 
in  the  stream  of  product  gases.  Some  particles  may  undergo  encapsulation  by  an 
oxide  crust,  followed  by  blow-out  of  the  magnesium  vapo;  (Fig.  la),  but  it  seems 


that  this  mechanism  may  prevail  only  at  the  walls  and  other  places  that  permit  the 
physical  entrapment  of  the  particles.  The  reaction  in  the  main  gas  stream  appears 
to  be  one  of  a continuous  boiloff  from  the  exposed  surface,  causing  each  particle  to 
trail  a cloud  of  vapor  or  oxide  {Fig.  lb). 

3.2.2  Determination  of  the  Magnesium  Droplet  Lifetime 


The  steady-state  diffusion  rate  of  the  magnesium  vapor  from  a droplet  through  a quasi  - 
stationary  layer  of  non-diffusing  gas  is  given  by 


dm  _ D 

dt  R T 

s 


4irr 


frs  ' 


Po> 


(15) 


where 


-jH-  = evaporation  rate,  mole/sec 


r = radius  of  fuel  droplet  at  time  t = 

Tg  - surface  temperature  of  fuel  droplet,  1363“  K 

p ~ partial  pressure  of  magnesium  at  surface  of  droplet,  1 atm 
s 

p = partial  pressure  of  magnesium  in  the  gas  stream,  0.12  atm  (see 
° Table  8) 

3 

p = density  of  magnesium,  1. 74  g/cm 


3/2 


D = 0.0043 


P Iv 


1/3 


+ V 


1/ay 


Vt 


M. 


(Ref.  12) 


(16) 


where 


and  = molecular  weights  of  air  and  magnesium 

V and  V.  = molecular  volume  for  M and  M, 
a b „ a b 

D = 358  cm  /sec 


Integrating  tbe  left  side  of  Eq.  {18)  from  to  0 and  the  right  side  from  t = 0 to 
t = t.  , *he  lifetime  of  the  diop,  one  obtains 


RT  r2  p 
s o s, 

2 D M (p  - e ) 

3 O 


0.5  x 10  sec 


(17) 


3,2.3  Calculation  of  Reaction  7one  Thickness 


During  the  time  increment  dt,  the  thickness  of  the  gas  layer  freshly  ejected  from  the 
soiid  surface  is 


(dx) 

e 


RT  m 
a 


vdt 


(18) 


where 


T = average  temperature  between  the  boiling  magnesium  (1363*K)  and  the 
a flame  temperature  (30004Kf  see  Table  3),  assumed  to  be  2000®  K 

m - moles  of  material  in  gas  phase  per  gram  of  solid  decomposed 
S =1.66  x lO-2 

= mole  f 'cMon  of  nitric  oxide  and  oxygen 
= density  /i  gaseous  mixture  (including  entrained  solids) 

v = regression  rate  (cm/sec) 
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At  time  t , the  thickness  of  this  oxidizing  gas  layer  is  reduced,  by  reaction  with  the 
magnesium,  to 


- (<*), 


(19) 


m , the  number  of  moles  of  magnesium  in  the  drop  at  time  t , is  given  by  integration 

of  Eq.  (15),  and  when  substituted  into  Eq.  (19)  gives  the  thickness  of  the  combustion 

zone , x : 

’ c 


x 


c 


2ab' 


5/2 


5 c 


(20) 


where 


a = 


R T m y v p v 
a gJN*g 

m 


c 


8 xD 

3 R T 
a 


Substituting,  one  obtains  x = 0.12  cm. 

V 

Equation  (20)  permits  the  determination  of  the  reaction  zone  thickness  as  function  of 
particle  size  as  shown  in  Figure  9.  The  finding  that  the  reaction  zone  thickness  is 
short  compared  with  the  depth  of  the  bullet  cavity  agrees  with  published  data  (Ref.  13) 
which  also  indicate  that  for  80  mesh  (175  p)  or  smaller  magnesium  particles,  the  burn- 
ing rate  is  not  affected  by  particle  size.  We  may  conclude,  therefore,  that  for  m x- 
tures  which  are  close  to  stoichiometric  proportions,  the  rate-determining  processes 
occur  in  the  solid  phase . 
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Section  4 
RESULTS 


Figure  10  shows  the  temperature  distribution  in  die  solid  phase  as  computed  with  the 
aid  of  the  computer  program  listed  in  AppK'Mix  A.  Figure  10  also  shows  that  the 
surface  temperature  of  the  so  ld  is  abovt  .e  melting  point  of  magnesium  because 
the  model  does  not  take  into  ac  ount  the  melting  of  the  magnesium  in  the  surface  layer. 
If  *’  s r fineraent  were  to  le  incorporated,  one  would  expect  the  surface  temperature 
to  lov'pr,  somewhere  b ween  1000  and  1300*  K.  Such  a lowered  surface  tempera- 
ture would  bring  "it  steeper  temperature  profile,  but  any  heat  contributed  by 
the  decomposing  PVC  would  tend  to  raise  the  temperature.  Because  of  this  foot, 
the  sicr  '1,r’cation  made  in  the  model  seems  to  be  reasonable. 

The  particle  size  of  the  magnesium  does  not  affect  the  temperature  profile  because 
the  mixture  was  considered  for  computational  purposes  to  be  homogeneous.  It  will  be 
noted  that  the  particle  size  of  the  magnesium  (0. 012  cm)  is  of  the  same  order  as  the 
solid  reaction  zone  in  the  vicinity  of  the  surface.  Because  the  melting  process  is 
fast  compared  with  the  decomposition  rate  of  the  strontium  nitrate,  one  may  indeed 
expect  the  surface  to  consist  of  molten  magnesium  particles  at  a temperature  which 
is  somewhat  below  the  boiling  point  { 1363°K).  This  finding  was  confirmed  by  a mo- 
tion picture  study  (Ref.  1). 

Figure  11  shows  the  computed  relationship  between  regression  rate  and  heat  flux  into 
the  solid.  The  model  does  not  account  for  the  spin  effect  on  regression  rate.  Further, 
the  effect  of  particle  size  (i.  e. , the  thickness  of  the  flame  zone  and  hence  the  bright- 
ness of  the  luminous  source)  is  contained  within  the  heat  flux  and  cannot  be  stated  ex- 
plicitly because  the  model  does  not  account  for  the  effect  of  the  reaction  zone 
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0 FUEL  RICH,  FLUX  250  CAL/CM2-SEC 
0 FUEL  RICH,  FLUX  500  CAL/CM2-SEC 
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® OXIDIZER  RICH,  FLUX  500  CAL/CM*-SEC 
0 STOICHIOMETRIC,  FLUX  250  CAI/CM2-SEC 
0 STOICHIOMETRIC.  FLUX  500  CAL/CM2-SEC 
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Fig.  10  Temperature  Distribution  in  Solid 
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Fig.  11  Computed  Reaction  Rates  as  Function  of  Heat  Flux.  Data  points  were 
' experimentally  obtained  at  zero  spin  rate 


thickness  on  the  radiative  heat  transfer  to  the  solid.  We  have  shown  (subsection  3.2), 
however,  that  the  path  length  of  the  fuel  particle  is  brief  with  respect  to  the  depth  of 
the  tracer  bullet  cavity  and,  therefore,  that  for  cases  of  particle  sizes  of  less  than 
20 flu  and  mixtures  not  very  f ar  removed  from  stoichiometric  proportions,  the  find- 
ings of  this  study  will  be  applicable. 

Table  13  lists  the  heat  fluxes  and  energy  partitions  for  die  cases  for  which  experi- 
mental burning  rates  were  available. 


Table  13 

CALCULATED  HEAT  FLUXES  AND  ENERGY  PARTITIONS 


Tracer  Mixture 

Regression 

Rate 

v,  cm/cec 

Heat  Flux, 
Fig.  11 
cal/cm2sec 

Energy 
Partition  See 
Appendix  B 

Oxidizer  Rich  at  n = 0 

28.8%  Mg,  71.2%  Sr(N03)2 

0.15 

125 

27% 

Stoichiometric  at  S2  = 0 
36.3%,  Mg,  63.7%  Sr(N03)2 

0.31 

230 

18.5% 

Fuel  Rich  at  f>=  0 
42.8%  Mg,  57.2%  Sr(N03)2 

0.41 

280 

20 

R-284  at  «=  40  000 

17%  PVC,  28%  Mg,  55%  Sr(N03)2 

0.27 

-220 

- 

If  the  model  selected  is  valid  in  the  presence  of  spin,  the  computed  results  shown  in 
Figure  11  can  also  be  used  to  study  the  effect  of  spin  on  the  burning  rate.  We  observe 
that  heat  flux  into  the  solid  increases  with  spin,  whereas  the  energy  partition  is  in- 
dependent of  spin.  There  are  two  reasons  why  the  flux  could  increase  with  spin:  (1) 
the  centrifugal  action  may  cause  the  ejected  particles  to  be  swept  across  the  surface 
and  so  reduce  the  thickness  of  the  boundary'  layer,  or  even  erode  the  shape  of  the 
solid  surface;  and  (2)  the  burning  mechanism  may  change  with  spin  as  by  the  breakup 
of  the  reactive  surface  leading  to  the  ejection  of  oxidizer  and  fuel  aggregates  into  the 
gas  stream.  It  may  be  helpful  for  the  understanding  of  the  burning  mechanism  to 


have  a motion  picture  study  of  the  burning  spinning  tracer  compositions,  although 
the  technical  problems  of  resolving  the  surface  phenomena  will  be  formidable.  A 
breakup  of  the  solid  would  give  rise  to  an  extended  reaction  zone.  It  may  be  possible 
to  assess  this  possibility  by  firing  tracer  bullets  without  spin  (as  from  a smooth  bore 
barrel  using  a sabot)  or  by  investigating  the  plume  generated  in  a laboratory  bullet 
spinner.  Such  studies  would  also  serve  to  confirm  the  findings  of  this  study. 
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Section  5 


RECOMMENDATIONS 


The  assumed  decomposition  model  of  the  solid  phase  might  bo  altered  to  allow  for 
some  oxidation  reactions  to  occur  in  the  solid  as  well  as  the  ejection  of  some  oxi- 
dizer aggregates  from  the  surface  prior  to  decomposition.  This  modification  would 
no  doubt  be  important  in  systems  that  contain  a liquid  phase  in  the  reaction  zone  such 
as  9odium  nitrate  (flare  composition)  or  in  which  the  gas  phase  reaction  is  of  second- 
ary importance  (smoke  generating  munitions).  A refinement  of  the  computer  pro- 
gram may  be  achieved  by  taking  into  account  the  melting  of  magnesium  at  the  surface. 
In  small -caliber  ammunition,  the  deviations  from  the  assumed  one-dimensional  ge- 
ometry may  become  significant,  and  corrections  for  the  heat  losses  to  the  container 
walls  should  be  incorporated. 

Application  of  the  model  to  other  types  of  munitions  that  may  operate  in  an  oxidizer- 
rich  regime  may  require  further  study  of  the  vaporization  law  used  in  this  study. 

The  presence  of  a diffusion  barrier  about  the  fuel  particles  (such  as  an  oxide  crust) 
will  then  become  a significant  factor  in  limiting  the  reaction  rate. 

The  mechanism  by  which  spin  induces  a change  of  burning  rate  may  require  further 
study.  If  the  spin  effect  is  due  to  a time-dependent  phenomenon  within  the  reaction 
chamber,  the  burning  rate  should  not  be  assumed  to  be  constant.  Except  for  the  ex- 
perimental study  required  for  an  improved  description  of  the  spin  effect,  all  other 
problems  are  amenable  by  an  approach  similar  to  the  one  chosen  for  this  study. 

A remaining  problem  area  is  the  modeling  of  the  energy'  partition.  This  task  would 
require  assigning  the  radiative,  convective  and  conductive  heat-transfer  components 
so  as  to  match  experimentally  determined  energy  partitions.  The  problem  is  made 
difficult  by  the  complex  radiation  phenomena  as  well  as  by  the  possible  deviation  from 
steady-state  conditions. 
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Section  6 

SUMMARY  AND  CONCLUSIONS 


The  modeling  of  tracer  reactions  permits  the  determination  of  the  parametric  influence 
of  composition  on  burning  rate  and  of  particle  size  on  the  reaction  zone  thickness . 
Neither  a rigorous  analysis  of  the  effect  of  component  and  mixture  characteristics  on 
trace  duration  and  trace  visibility  nor  the  mechanism  of  the  effect  of  bullet  spin  on 
burning  rate  was  included  in  this  study.  The  tracer  reaction  appears  to  be  completed 
within  the  bullet  cavity  so  that  the  composition  of  the  exit  gases  is  accessible  from 
equilibrium  calculations . The  influence  of  tracer  compositions  on  the  composition  of 
exit  gases  was  determined. 

Thermal  transport  into  solid  tracer  mixtures  was  found  to  be  relatively  insensitive  to 
changes  in  composition.  It  was  shown  that  all  thermal  processes  that  occur  in  the  solid, 
such  as  dissociation  and  melting,  occur  with  200  p from  the  surface,  i.e. , within  a 
depth  of  no  more  than  two  fuel  particles.  The  burning  rate  varies  approximately 
linearly  with  heat  flux,  increasing  with  increased  fraction  of  fuel.  At  the  extremes  in 
composition , this  generalization  is  not  expected  to  apply. 

The  model  at  present  cannot  predict  the  burning  rate  a priori  without  an  experimental 
basis  for  the  calculation  of  the  energy  partition  which  must  account  for  the  effect  of 
spin  on  the  reaction  mechanism  and  for  the  relative  importance  of  heat  transfer  modes 
within  the  bullet.  At  zero  spin  it  was  found  that  the  burning  rate  is,  tor  fuel-rich  sys- 
tems, Insensitive  to  the  fuel  particle  size,  and  that  this  finding  is  in  agreement  with 
previous  studies.  The  methods  employed  in  this  study  are  expected  to  be  of  value  in 
the  study  of  other  pyrotechnic  devices  such  as  fumer  bullets,  flares,  and  smoke 
canisters. 
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Appendix  A 

COMPUTER  PROGRAM  FOR  THERMAL  TRANSPORT  IN  THE  SOLIDS 


A,  1 ORGANIZATION 


After  reading  in  the  necessary  input  data,  the  program  sets  up  the  dimensionless 
variables  and  parameters,  through  the  use  of  dimensionless  groups. 

Next,  the  input  data  and  dimensionless  groups.  are  written  out. 

The  calculation  is  then  carried  for  successive  time  increments.  At  the  end  of  each 
time  increment,  various  dimensionless  quantities  are  displayed  for  each  distance 
increment. 

A.  2 INPUT 


C 

D 

K 

Q 

Z 

E 

HD 

TOD 

DXD 

P 


= heat  capacity 

- density 

- beat  conductivity 

= heat  of  endothermic  reaction 
= pre  exponential  factor 
= activation  energy 
= heat  flux 

- initial  temperature 
= distance  increment 

= ratio  of  time  increment  to  distance  increment 
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EPS 

NL 

MXT 

KU 

ITERMX 


small  number  indicating  when  decomposition  is  complete 
number  c.“  distance  increments 
maximum  number  of  time  increments  allowed 
geometry  factor  (zero  for  one-dimension  geometry) 
maximum  number  of  iterations  allowed  in  the  calculations 


A.  3 OUTPUT 


For  each  time  increments  or  a multiple  number  of  time  increments,  the  program 
prints  out  the  following  dimensionless  quantities  for  each  distance  increment  or  for 
predetermined  distance  increments: 


J 

= 

time  index 

TIME 

= 

dimensionless  time 

I 

distance  index 

Dist 

= 

dimensionless  distance 

T 

» 

dimensionless  temperature 

F 

= 

fraction  of  nitrate  undecomposed 

When  decomposition  is  complete  in  a distance  increment,  a message  is  printed  out  to 
this  effect.  The  output  table  is  then  presented  with  that  increment  omitted. 

Dimensionless  time,  distance  and  temperature  can  be  converted  to  their  dimensional 
counterparts,  using  multiplication  factors  given  in  the  printout. 

Regression  rate  can  be  obtained  simply  by  dividing  a burnec’-through  distance  by  the 
corresponding  time. 
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1* 

2» 

J» 
4* 
!>• 
6 • 
'/• 
8* 
V. 
10* 
u* 
« • 
13* 

14» 

IS* 

u» 

17. 

18* 

tv* 

20. 

21* 

22* 


2S« 
26* 
27* 
28* 
29. 
30* 
31  • 
3J* 
33* 
34. 
3S* 
36* 
3/4 
38* 
3V* 
40* 
41* 
42* 
43* 
44* 
4S* 
46. 

47* 

48* 

49* 

30* 
51* 
S2* 
53* 
54. 
5S* 
!»• 
57. 
58* 
5 V* 
*0* 
it. 
62* 
63* 
64* 
6S« 
664 
67* 


C 


c 


c 


c 


RC*L  X - • 

DIMENSION  Ti20U,4>,  P«2O0‘,4> 

DIMENSION  AlFAiZOM*  BETa<2o0>*  C*M4t2Cn).  PHt!200>,A!200)iPr**20Q> 

DIMENSION  Dfl.T*«200»,  Ct 200 ) 
too  P'CH  lATlflF  10*  4 > 

102  FEH-<4T(B|10) 

2UU  *C»*'»TtlHl*  4 ho  • i E 1 0 1 4 5X»  4H0  « l El3.4t.5Xt  4HK  * , El0*4,. 
*Sy,  4 HO  * . tlU,4.  5x»  4H?  . , Eir.4|  Sx»  4H£  » l Et0*4) 

203  FOR'ATUHO.  SH»10  • t r.10,4,  5X,  6HT00  * » Eto*4*  5x*  6M0X0  • • 
•fcio.4*  5X.  * HP  * . *U!3,4,  Sx,  6HEP5  * * EtO.«> 

204  FOHM*T(1hO,  SHNL  ■ , 1 4 , <1X,  6HMXT  ■ i 1 4, . 5X* . SHNU  • .MA^-SX.- 
*»H1T£KHX  « , 14) 

206  FCH’UT<iH6.  9HIEMFCT  • « F10,4,  SX.  •MTIP^T  • » E10.4,  Sx, 
•VlOSTfcT  • , HO, 4,  SX,  4wB  • , E10,4»  3X*  4«H  * , E10.4) 

2r8  FCH  UT(//«HP|  >X,  7MJ  l^DEX,  4X.  THIt.  TIHEW  3X,  7P1  |NDEx.  4v, 
•7|iUL  CIST,  7X,  7H0L  TEMP,  7X,  10HDL  UNREAL) 

2lU  FoH'ATUhO,  Sx.  J3,  5X,  £<0,4,  5X.  (3,  9X.  3(El0i*.S*h 
212  FON'iATt  1H0,  1(1X,  7H1TER  » * 13.  10X*  "HRAtHX  • • E10,4) 

214  FCHUT<///iHo,10x,21HH£w  L»rEH  IS  DEPLETE Jt) 

.....Input  data....* 


KHOlS.lon)  t*  0,  K,  0|  z‘,  E 
HE*rilS,10O>  H(),  TOO*  0X0,  p,  EPS 
HMi)i5*lg2)  NL.  H*T,  HU*  fTERPX 
WBlTt<6, 20(11  C,  0,  X,  0.  t,  E 
W«lTE<6.2n2)  HO*  TOO,  0*0*  P*  EPS 
WRJT£«6,204)  NO,  HXT,  NU,  ]TfRMX 
•••••|N|TMU24T10H*9**« 

H « 2*0 
TEHFCT  , R/F 
TlMfCT  « R«r)«Z/(VF. 

OFTrc'  * SCRT<K.o*a»Z/K7£) 

B « H*0/c/E 

H * H0.SvJRT(H/U/Q/2/K/E1 
TO  * TfcMF CT. 1 00 


OX  • DSTFn*D*U 
ot  « p*nx**? 

Pn  ■ 1.0  ♦ P 
PM  • 1.0  - P 
PC2  f P/2,0 
PpXu2  « 2,0»P*UX»H 
NIM1  * lit  - 1 

WHHE(6,206)  TcMfCT , TJUp^T,  DSTFcT,  B,  H 
UC  300  1 « 1,  NL 
Tt 1.1)  * TO 
F (1 , 1 ) * 0 
300  CCNTJNUC 
T l Hi;  « 0,0 
IS  ■ 1 

.....  START  00  LOOP  FOR  TIME  INDEX  U 4.*.# 

Do  ^32  J * 1,  MXT 
JLHITt  « 0 

IMJ  ,LE , 5 ■ OH , HOOIJ.S)  ,EQ,  0)  jRHJTE  * 1 
IF( JWR  1T£  ,ru.  1)  WN JTE( 6,206) 

Tint  * TJHt  ♦ UT 
N » M00(J-;,4)  * 1 
M » M0UIJ.4)  ♦ 1 
NN  » H00(J*1*4)  . 1 
HF  * HOO(J+2|4)  ♦ 1 

*4*»«  F1HST  STfcP  IN  APPROXIMATE  SOLUTION  **••• 

DO  304  1 * IS,  NL”i 

*!l>  » F!.OaT(NUI/?,0/FIO*tC!»1) 

pptll  • exf(«i*o/t<i,nj»/r.dt 

IF<  I iCT,  IS)  go  TO  }C6 
Bet*i i ) • ppi 

GAHA( I ) * .P 

PhH|>  • PMl*TlI,N|*p.Til*i,N)*POXH2.Pll,4).pP<lV 
DELTA! |)  • CAPA! 1)/*eTA< I) 
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! 

t 


IH* 

6V* 

TO* 

71* 

7 a* 

73. 
71* 
72>* 
7ft* 
77. 
7H* 
7 V* 
BU* 
Bl 
82* 
83* 
#4* 
6b* 
86* 
87. 
8H. 
8*#. 
90* 
91* 
92* 
93* 
94. 
9b. 
9b* 
97* 
9H* 
9 V* 
100* 
101* 
102* 
103* 
104* 
10b* 
10b* 
10?* 
108* 
10** 
110* 
111* 
il'<?* 
113* 
114* 
lib* 
lib* 
117* 
ill* 
11’* 
120. 
121* 
122* 
123* 
124* 
12b. 


cm  * PHI  { ? )/b£T At  ? > 

Gr  TO  304  ... 

3t<8  AileC!)  * Pf>2*tl,n-M  1>> 

BfIM!)  * “PHI 
ca  .*m  * PC2*ii,n*»<t)S 

Hwicp  * -pn2.ii, o.Mm.Ttj.i.Ni-pHi.rn.Ni.Poa^a.o.Ad n* 

*T< I+liN)  + F(!.N).PP(!) 

DlLTAtl)  « I )/(yFT*C  t I l.QELTAt  J.lJ  > 

G<  1 1 = (PHim*ALF*(I).Gu.1))/(BEY*m*A.F*(!>*DEUTAtl«in 

3c-4  continue 

ALF4|NU)  * ppi 
arT/jNu  ■ -h 

PRKNL)  * Pf’l*T(NL.N>4P*T<NL-l#Nl-F<NL,N)*PPjNL> 

G(Ni>  = <PHlcNL)-*LFA^Ui.C{NL.Uw<BETA(9L)*ALFAlNL>.0EuU+Nt.a*->- 
C GETTING  IKE  *PpROXf*ATE  SOUJTtON 

00  306  | * NL,  IS.  -1 
If  1 1 ,LT i NU  GO  TO  310 
Tt|,M>  « G<  1 > 

1 r f tC I * m)  .lt«  To)  T< 1*n»  • TO 
Flt.H)  « Ft  l|N)*U,0«PPt  |>) 

GC  TO  306 

310  Ti  1 = G(l>  * OULT*(  I )*Tt  I*1»M 

(KTIIfH)  ,LT.  TO)  Tll.tl)  « TO 
FU.rD  * Ft  !,M)*tl,0«PPl  1)1 
iptrl I »m*  ,LT,  0,0)  FU»M?  « 0,0 

306  CONTINUE  . 

IMJ^RITE  , NEi  1)  GO  TO  318 

C WHSTIMG  OUT  THg  APPROXIMATE  S0UUT1ONI  *•••• 

00  312  I ■ 1,  NU 
D 1 5T  = <I-l>*OX 

IF 1 1 ,LT , IS)  wRlTfct6»210)  J.  T , DIST 

|rt  I »GE,  IS)  WRlTtlft.ZlO)  Ji  - I*  DISJ,  TU#Ml,.jrl4,R) 

312  CCNTJNUe 
33 « CONTINUE 

C »••••  FIRST  STfcP  JN  IMPROVING  THE  UTJON.  *•••• 

UC  320  I « ts,  N'u 
TCl.NN)  « T< | i M) 

320  CONTINUE  

C .....  bUGIN  THE  ITERATION  ••••• 

JTtR  * 0 

328  ITEM  * ITER  ♦ 1 

Nc*r,  * o 

K*T  >iX  * 0.0 

UC  314  I « IS,  NL  -.  . 

PX  * DT/B*tXP(-2,0/(T< I»NN)*T«  l,N>>> 
lr< I , GT , IS)  GO  TO  316 

T ( I « rift ) * l,07PPl*(PMl*Tt  f.N)'P*(Tl  I*1.NN)*T(  I*1iN))*PdXh2 
•*il,0“tXP(-PX)).Ft  1 1 N> ) 

GO  TO  322 

31b  IF  • I ,EJ.  NU?  GO  TO  J18  _ _ 

Tt  l, MM)  ■ 1,0/PP1*(PM1.T(  i,*i)4P02.(U,0*A(  i>  >•<  Ti  J*1,NN)*T<  1*1, N)j 
*♦(  1. Q’A(  I ))*<Tt  l-i.NN)*Tt  l-i,Nl>)«(i,0-EX:>lt»PX))*F(I,N>) 

GO  TO  322 

318  T(J.MM)  * l,0.'PPl«(PMl*TI|,N)*P.(Tt|*l,NN)*T(|"l»N)>« 

•<1.< -EXP(-PX) )»F( I . MJ ) 

c •«**•  1 STATEMENT  for  testing  »*•* •••.•••••*•••••**. 

T(l.MM)  t TO 

322  RaT  ■ (T!l,MM)-T(IrNN>>/T!!,NN) 
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126* 
12  f* 
12H* 
129* 
130* 
til* 
132* 
133* 
134* 
13b* 

136* 
13 /« 
13«* 
i 3M* 
140. 

141* 

142* 

143* 


144* 

14b. 

l4h* 

14/* 


14H* 

149* 

ibi]* 

lbl* 

lb2* 

lb3* 


154* 
lbb* 
1 b'  » 
lb  - 
lba* 
lbv* 
160* 
161* 


162* 


ABS.'AT  = AaS(KAT) 

IfUgSHAl  , GT , EPS)  nCAG  = 1 
If  ( 'vSSKaT  ,G3,  RaTMV)  Ri T*X  s AijS^AT 
pyx  S UT/e.EXPi-i.O/fTd^Hj  + TU.N))) 
ft  l, MM)  a Ft  l,N)*F.XP(-PXx) 
irii-(I'MM)  , L,T » 0,0)  FH.MH)  a o,n 
3i4  cr^Ti^Uf. 

C *****  IMIS  COMPLETES  CNE  fTtHATlOr.  ***** 
if.L.  » 0 
JWK  < * f) 

Ur  >26  I s lb,  \L 
T t I » Ml-1 ) a T ( 1 , hM ) 

KAl  jr  = y ( i , rft)/5 

iFt'.-ATf  ,LT , i,r,t*o3)  Is0  s ISO  ♦ , 

IF(:,aTF  ,LT,  l,tit-U3j  jWp*  = 1 
326  Cpi\T  1 91ir 

IF  f I T£K  ,LE*  IIEHmX  NOaG  ,E0.  U GO  Tc  32fl 

IrtjWNiTt  ,N't»  1 .AMq,  JwqX  ,Nfc.  ))  GC  TO  34o 
IF(  JyRITfc  .ft),  j , A 10 » -July  , Nt*  «)  GG  TC  ip, 
WF.iTE(6,2i4)  * X 

l **•»*  WBlTjr^u  OUT  Tng  IMPROVED  -Ol_UTjC\  ***** 

3C4  WMict6,2i2)  HErs,  TaTmX 
DC  >30  I a l,  iml 
UIS1  * ( i»i)*nx 
I r < I • LT , IS)  00  10  334 

•Jf  ite<6,2iu>  J»  time,  j.  pist.  fIi.m) 

GC  T0  330 

334  WF  1 te.<  0,210)  J»  TIME,  I» 

330  Cfi'TIWUp 

o *****  STaPT  THt  NEXT  T I Mj:  CYCLE  ***** 

340  If  = lb  * I CD 

UC  .542  I a 1,  NL 
T ( I , H ) = T ( I » HM  > 

Ft  l . D * Ft  I, MM) 

34  2 COTJNUF 
332  CC'iTI  JUL 
tM-1 
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Appendix  B 

CALCULATION  OF  ENERGY  PARTITION 


Heat  Generation,  Stoichiometric  Case 

5 Mg  + Sr  (N03)9  5 MgO  + SrO+  N2 

AH  =■  1880  cal/g 
R 

Rate  of  Heat  Generation,  Case  of  Zero  Spin  (see  Tables  12  and  13) 

cal 

AH  D x p x v = 1235  — y 

R r sec  cm* 


Heat  Flux  into  Solid  (see  Fig.  11) 


cal 

H — 230  2 

sec  cm^ 

Energy  Partition:  230/1235  = 18. 5% 


Similarly,  for  Spin  Rate  of  40, 000  rpm 

H = 340 

340 

Energy  Partition:  TZ7T  = 18.5% 
1800 


It  must  be  remembered  that  Fig.  11  was  not  based  on  a model  which  included  a spin 
effect.  It  is  possible  tnat  a refined  model  would  not  yield  the  same  conclusion,  i.e. , 
that  the  energy  partition  is  independent  of  spin. 
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